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Abstract 

In a nanopillar with dipolarly coupled vortices, we present an experimental and simulation study 
to understand how the interplay between the bias field and spin transfer torque impacts reversal 
of the vortex cores. We find that, depending on the current values, vortex cores might experience 
different physical mechanisms for their reversal, namely a static or a dynamic switching. We believe 
that our results might be useful in the context of vortex based non volatile memories, as a current 
controlled selective core switching is proposed . 
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Remarkable static and dynamical properties of magnetic vortices, and particularly the one 
of its cores jl, 2] has recently motivated numerous experimental and/or theoretical studies 
to investigate the behavior of a single vortex under the influence of current induced torques. 
Notably spin torques associated to rf currents have been proved to be an efficient way to 
reverse a vortex core. As for the vortex dynamics, very large amplitude of gyrotropic motion 
of vortex core can be obtained due to spin transfer torques, making vortex basic systems 
very interesting for the investigation of fundamental aspects of vortex dynamics 0, 0] and 
even chaotic behaviors [5] but also for their implementation as highly integrated local rf 
oscillators in telecommunication systems jfj]. 

One of the most interesting properties of magnetic vortices are their high stability in 
geometrically confined systems, explaining why these magnetic objects are good candidates 
as a new media for non- volatile magnetic memories in which information is stored either in 
their core polarity or chirality (3, S|. Thus it is of major interest to identify the physical 
mechanisms, often related to the static and dynamical magnetic properties of the material, 
that are controlling the core reversal processes. These issues are primordial for the optimiza- 
tion of the writing procedure but also for increasing the final speed of these vortex based 

n n 

memories. First studies focused on static field induced core switching [2J, |9j , in which the 
vortex core breaks down and form back with opposite polarity. This mechanism leads to 
relatively large switching fields that are mostly material dependent (and thus not easy to 



modify). Recent 
or an rf current 



14 



resonant core reversal under the presence of either an rf field 1CH13| 



15] have been proposed to lower their switching fields but also in the 



perspective of memory devices, to provide a selective tool to reverse a single core. Finally, 
core reversal have been predicted using injection of dc current and the spin transfer torque 
thereof 

Our main objective is to investigate how, in a system with two coupled vortices, dif- 
ferent core reversal processes can coexist, namely static reversal at low dc currents and 
dynamic reversal at large currents, and more importantly how the current induced dy- 
namics of the coupled vortices can allow to switch between the two processes. In our 
system i.e. a spin valve nanopillar with one vortex in each of the two magnetic layers, 
we have performed a detailed experimental study supported by micromagnetic simula- 
tions to understand the evolution of the core switching fields under the action of a large 
dc current (i^) and a perpendicular applied magnetic field (H perp ). The sample stack is 
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FIG. 1: Resistance versus H perp at Idc = 11mA measured on 0l2Onm nanopillar. Vertical dashed 
lines indicate switching events, and vortices polarities are sketched for each region: top (red) 
correspond to thin layer's and bottom (blue) to thick layer's. 

Cu(60nm)/Py(15nm)/Cu(10nm)/Py(4nm)/Au(25nm) where Py=NigoFe2o, in which circu- 
lar nanopillars of two different diameters, 1 20nm or 200nm, have been patterned by standard 
e-beam and ion etching lithography process. Such spin transfer oscillators based on coupled 
vortex dynamics have demonstrated a record linewidth in the tens of kHz range even at zero 



field, making them very promising for the real development of spin transfer rf devices [191120]. 

Our experimental procedure is to sweep the perpendicular applied field H perp from large 
negative to large positive value at different injected dc currents Idc- In the starting configu- 
ration, the two vortices have their core polarities oriented parallel (P) and negative, and with 
identical chiralities (imposed by the Oersted field symmetry). This experimental procedure 
was applied for both 0l2Onm and 02OOnm circular pillars. We also performed micromag- 
netic simulations 33| on the 02OOnm nanopillar. The simulations take into account the fact 



that both magnetic layers in the stack can act simultaneously as a spin polarizer and as a 
free layer when considering spin torque action. The current flow is assumed to be uniform 
through the pillar and the Oersted field is taken into account. A constant small in-plane 
field H x =50 Oe is introduced to break the system's symmetry, but is not impacting the 
gyrotropic frequencies or any other dynamical properties. For each field and current value, 
we calculated the vortex dynamics of the nanopillar, with the initial state corresponding to 
the quiescent state at this field (i.e. with no current). 

In Fig. [U we show a typical R(H perp ) curve (measured on the 120nm diameter sample at 
Jrf c =llmA), that describes how the core switching fields are precisely detected. At H perp =- 
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4kOe, the two vortex polarities are parallel and oriented downwards in our field convention 
(Pthin = Pthick — — 1)- Then we sweep the field towards positive values. In this field range, 
a smooth parabolic variation of R is observed that is not related to a change of the vortex 
core configuration but rather attributed to the reversible variation with H perp of the out-of- 
plane magnetization component in the two vortex tails. The first direct evidence of a core 
switching event in Fig. [TJis observed for H perp =2.6kOe at which a sharp resistance variation 



occurs. We recently proved 



191 ] that this resistance change comes from the fact that a 



large amplitude spin transfer gyrotropic mode is generated when the two vortex polarities 
are AP, increasing the core-core distance and subsequently the mean resistance. In such 
configuration with AP polarities, a large rf peak can be detected with a frequency in the MHz 
range corresponding to the gyrotropic mode of the coupled vortices, typically /=l.lGHz at 
7<fc=llmA for 0l2Onm, and /=720MHz at Jd c =20mA for 02OOnm. An important feature 
of this vortex dynamics is that the main characteristics of the excited coupled mode are 
close to the one of the isolated vortex in the thick layer. Consequently, the gyration of thin 
layer's vortex core will be neglected in a first approximation. 

We emphasize that two sources of spin polarization can promote large amplitude spin 



transfer excitations 



211 ] of the vortex core in the thick layer. A first one is due to the 



perpendicular spin polarization associated to the thin layer's vortex magnetization t 
develops a perpendicular component under the application of the perpendicular bias field 



hat 
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221 ] . The second one is associated to the in-plane circular magnetization distribution of the 



thin layer's vortex that can equally contribute to the vortex excitation [23:]. A precise 
description of the evolution of these two contributions with bias field is out of the scope of 
this paper. However we can notice that they can add up and thus increase the gyration 
amplitude of the core when the magnetic field increases in the direction of the thick layer's 
core. A mean to identify which core has reversed is to measure the slope of the linear 
evolution of the frequency with H perp , HH perp ), because it is directly related to its core 
polarity sign : Pthick = sign(df / dH perp ) [24j. This feature allows us to ensure that the first 
resistance jump in FigJTJis related to the switching of the vortex core in the thin layer. In 
the field range where the two polarities are AP, the resistance remains on this plateau with 
large values, until it decreases sharply back to the low resistance level at H perp =3.1kOe, at 
which the core direction of the vortex in the thick Py layer eventually reverses. Note that 
each core switching is an irreversible process. We have repeated such R(H perp ) curve for Id c 
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going from 10 to 17mA. In addition, we have performed similar measurements for a 0'2OOnm 
nanopillar with going from 14 to 40mA. For each curve, the values of the two switching 
fields are extracted and linked either to the thin or the thick layer's vortex core switch. 

In Fig |5^a) and (b), we display the evolution of the two switching fields as a function 
of Idc for each pillar diameter, showing very similar trends for both nanopillars. Indeed, 
the variation of the two switching fields, H™ in and H™ ick with !<& can be divided in two 
regimes. First, at "low" currents {Idc < 14mA for 0l2Onm and Idc < 23mA for j2S"200nm), 
the thin layer's vortex core is found to switch first for both diameters. Furthermore, the 
switching field H™ in is almost constant on the whole low current range, whereas Hf^ ick 
decreases significantly to reach the values of H™ in at the end of this first region. As we will 
discuss later, these features are related to static core switching. In the second regime, at 
"high" currents {Idc > 16mA for j2fl20nm and Idc > 33mA for £f200nm), we find an opposite 
behavior, as the switching field of the thick layer's vortex core is always found to be smaller 
than the one of the thin layer's vortex. This is indeed indicative of the influence of the spin 
transfer induced vortex dynamics on the core reversal process. Moreover, H™ ick decreases 
strongly with J^ c , passing through zero-field value and reaching negative field values for 
Idc > 37mA in the nanopillar with £f200nm (see Fig^b)). This is a direct indication of the 
dynamic switching of the vortex core, as we discuss in more details below. On the contrary, 
we notice that H^ in remains constant in the whole large current range but with a mean value 
strongly reduced compared to the one in "low" currents. Between these two current regimes, 
we find a range where no change of resistance in R(H) curves is detected, indicating that 
both vortex are switching simultaneously at the timescale of the measurement. Indeed, the 
two switching fields are probably close enough so that any switch of one of them immediately 
triggers the switching of the second one. 

In Fig. MJd), we present the simulation results for the evolution of the two core switching 
fields as a function of Idc- Depending on the field and current combinations, either static or 
dynamic switching of the thick layer's vortex core might exist whereas only static switching 
of the thin layer's core is predicted (it is explained below how to differentiate between static 
and dynamic switching events in the simulations). We find in these simulations that both 
cores experience a static switching at low I dc , and the dependence of the switching fields 
for both vortices on Id c is very weak. At large Idc-, we find the static core switching of the 
vortex in the thin layer and the dynamic core switching of the vortex in the thick layer. For 
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FIG. 2: (a) and (b) : Experimental switching fields of the thin layer's (solid black circles) and the 
thick layer's (solid red squares) vortex cores (a) for the 0l2Onm pillar, and (b) for the 02OOnm 
pillar, as a function of (b) : Switching fields of the vortex cores for the 02OOnm pillar 

obtained from micromagnetic simulations (open symbols), (c) Corresponding simulated velocities 
of the vortex cores prior to the switching events for the thin (black circles) and thick (red squares) 
layers'. 
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this latter case, the switching field Hf^ ick decreases strongly with I^c and it even becomes 
negative. All these observations are in remarkable agreement with the experimental data. 
In the simulations, the two vortices' cores are always found to switch independently (no 
"middle" current range), however this may be related to the absence of structural defects 
and grains in the simulated systems and to the fact that the simulations are done at T = 
OK. In Fig[2[c) , we show the velocities of the cores prior to their reversal extracted from the 
simulations [34| . 

Recently, several studies based on micromagnetic simulations emphasized important dif- 
erences between the static and dynamic switching processes. For example, Thiaville et al. 



25j have shown that the "static" core reversal is mediated by a formation and propagation 



of a Bloch point (magnetic monopole) within the center of the vortex core. In the continu- 
ous micromagnetic limit, the switching field, which is the field at which the Bloch point is 
introduced into the core, is infinite. There is no precise numerical solution of this problem; 
the apparent field of the static core reversal yielded by the simulations is mesh-dependent, 
thus has a vague physical meaning. However, as was also shown in the Ref. 25| . the energy 
barrier for the core reversal reduces with an external applied field. This implies that the ther- 
mal fluctuations will result in the static core reversal when the barrier becomes comparable 
with the thermal energy (/cgT). The second process, that we referred to as "dynamical" core 
reversal, takes place through the deformation of the vortex shape because of the onset of a 
dynamical gyrofield oriented oppositely to the initial core polarity. Consequently, a vortex- 
antivortex pair is nucleated leading to the annihilation of the ori gina l vortex and eventually 



the presence of a single vortex with opposite core polarity [13 
this process (formation of the vortex-antivortex pair) can be correctly accountec 
simulations, moreover leading to an accurate prediction of the critical velocity 



13, 120, 123| ■ The onset of 
for by the 



Note 



that the critical velocity for the dynamical core reversal is independent on the size or shape 
of the sample j^, but only on the material parameters and perpendicular bias field. At 



zero perpendicular bias field, the critical velocity is given by [30J]: v cr i ~ I.667V/I, valued 



Dias 



3l|. 



at 306m/s and 320m/s for the thin and thick layer's respectively. When the magnetic 
field H perp is applied, the critical velocity changes approximately linearly with the field 

These facts provide us with a deeper insight into relationship between simulations and 
experimental data. In the "low" currents range, the spin transfer induced dynamics are 
associated with small orbit radius, and hence small velocity. As a consequence, the condition 
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that the vortex core of the thick layer reaches its critical velocity is never fulfilled for any 
current-field combination. The two reversals then correspond to core breakdown, starting in 
the core center without formation of a vortex-antivortex pair. Indeed, a potential dynamic 
reversal at the same field would require much larger velocities that the actual ones reached 
by the cores prior to their static switching (see Fig 12(c)). [35]. Moreover, we find that, 
at any field value, the core profile is narrower for the thin layer's vortex. Consequently, 
it results that the energy barrier for the static reversal is smaller for this layer, making it 
easier to switch by thermal fluctuations 25] . In simulations, the static switching fields are 
independent of !<&.. We believe that the experimental decrease of H^ ick with (see Fig. 
[2(a) and (b)) is dueto its finite velocity before the switching, that leads to a narrowing of 
the vortex core. This feature reduces the energy barrier for the switching while the thermal 
fluctuations increase because of the heating effect by the current. 

For "high" currents, while sweeping the field the system reaches a point at which the 
evolution of the thick layer's core velocity crosses the critical velocity for a given H perp 
before the vortex core breakdown field Hf^ iC is reached. Therefore, its core eventually 
switches before the thin layer's one. As shown in Fig [2(c), the reversals of the thick layer's 
core then happens at velocities that strongly increase when 1& C increases, corresponding to 
the decrease of the observed switching field with Id c , in agreement with predictions by A.V. 
Khvalkovskiy et al (3l|. Interestingly, this core dynamical switching can happen even for 
negative field values, i.e. an applied field oriented in the direction of the initial core polarity. 
This is then an experimental confirmation that spin transfer core switching induced solely 
by a dc current is possible even at zero field. In the high-current regime, when the field is 
increased further, the "static" critical switching field of the thin layer's core is eventually 
reached. Similarly to the "low" current regime, the switching field of the thin layer's vortex 
due this "static" process is almost independent of I<i c . Notably, we find different values 
for the static switching field of the thin layer's vortex between the two current regimes in 
the experiments. This large difference might be related to the dipolar interaction between 
the two cores: in the "low" current region, the thin layer's core switches being in the (P) 
configuration, whereas in the "high" current region, the thin layer's core switches being in 
the (AP) configuration. In conclusion, we present a combined experimental and simulation 
study of the variation of the switching fields of vortex cores as a function of an injected dc 
current in a spin valve system containing two vortices, one in each magnetic layer having 



S 



different thicknesses. We identify the mechanisms responsible for the core reversal under the 
combined action of spin transfer torque and perpendicular field opposite to the initial core 
polarity. At low currents, both cores reversals occur through a static reversal of the vortex 
core. In this regime, the switching fields are defined by the intrinsic material parameters and 
their dependence with 1^ is weak.. At large current, the large amplitude gyrotropic motion 
induced by spin transfer allows for the thick layer's vortex to reach the critical velocity 
for a dynamical reversal, mediated by the creation of a pair of vortex-antivortex. In this 
regime the switching field of the vortex core has a significant dependence on the current. 
Our results are thus interesting in the perspective of development of non volatile vortex 
based memories as we demonstrate the possibility to obtain a selective switching event in a 
two-vortices based devices, simply through the application of a combination of dc current 
and pulsed field. 
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ll-NANO-0016), EU grant MASTER No. NMP-FP7-212257, RFBR grant 10-02-01162 is 
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